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ABSTRACT 

This paper presents late-time near-infrared and Spitzer mid-infrared photometric and spectroscopic 
observations of warm dust in the Type Iln SN 2005ip in NGG 2906. The spectra show evidence 
for two dust components with different temperatures. Spanning the peak of the thermal emission, 
these observations provide strong constraints on the dust mass, temperature, and luminosity, which 
serve as critical diagnostics for disentangling the origin and heating mechanism of each component. 
The results suggest the warmer dust has a mass of '^ 5 x 10"** M©, originates from newly formed 
dust in the ejecta, or possibly the cool, dense shell, and is continuously heated by the circumstellar 
interaction. By contrast, the cooler component likely originates from a circumstellar shock echo 
that forms from the heating of a large, pre-existing dust shell ~ 0.01 — 0.05 M© by the late-time 
circumstellar interaction. The progenitor wind velocity derived from the blue edge of the He 1 1.083 
fiia P Gygni profile indicates a progenitor eruption likely formed this dust shell ^100 years prior to 
the supernova explosion, which is consistent with a Luminous Blue Variable (LBV) progenitor star. 



Subject headings: circumstellar matter 
— dust, extinction — 



1. INTRODUCTION 



— supernovae: general 
infrared: stars 



supernovae: individual: SN 2005ip 



SN 2005ip was discover ed in NGG 2906 ( d « 30 Mpc) 
on UT 2005 November 5 (|Boles et al.ll2005f ). Early opti- 
cal spectra suggested th e discovery occurred a few weeks 
following the explosion (jModiaz et al.ll2005l ). The super- 
nova is a Type Iln event given the deve l opment of narrow 
Ha emission lines (jSmith et al.ll2009bD . iFox et all (|20090 
discovered a late-time near-infrared excess that has per- 
sisted at least two and a half years post-discovery. 

Late-time infrared emission typically indicates the 
presence of warm dust. The origin and heating mecha- 
nism of the dust, however, is not always well constrained. 
The dust may be newly formed or may have existed at 
the time of the supernova. If newly formed, the dust may 
condense from refractory elements in the expanding su- 
pernova ejecta or in the cool, dense shell of post-shocked 
circumstellar gas lying in between the forward and re- 
verse shocks (like the post-shocked environment trail- 
stellar wind collisions in WR binary systems (jUsovl 
In both cases, several heating mechanisms are 
possible, including radioactivity, optical emission from 
circumstellar interaction, and collisional heating by hot 
gas in the reverse shock. 

Alternatively, pre-existing dust may be collisionally 
heated by hot, shocked gas or radiatively heated by cither 
the peak supernova luminosity or late-time optical emis- 
sion from circumstellar interaction. In the latter case, an 
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'IR echo' is evident due to light travel time effects. Multi- 
ple scenarios can contribute to the late-time infr ared flux, 
as in the cases of SNe 2004et (iKotak et al.|[200 9. Suger- 
man et al. in prep) and 2006ic (iMattila et al.ll2008ll . 

Type Iln supernovae, na med for their "narrow" emis- 
sion lines (jSchlegell Il990( ). are more often observed 
to exhibit late time infrared emission associated with 
warm dust tha n any other superno va su bclass (e.g 
Pastorello et al.' '2002': 'Gcrardv et al.' '2002'; 'Pozzo et al 
2004; Fox et al. 2009; Smith ct al. 2010; Miller ct al 
2010allb l. Disentangling the origin and heating mech- 
anism of the dust can yield important diagnostics con- 
cerning the circumstellar medium, supernova progenitor, 
and explosion dynamics. For example, if the dust is 
shock heated , the o bserv ed dust temp erature yields the 
gas density ()Dweklll987t iDwek et al.ll2008[ ). which can 
be used to trace th e progenitor's mass loss history (e.g., 
iSmith et al.ll2009bl ) . If the dust is newly formed, the ob- 
served dust mass can be compared to models that predict 
supernovae as primary sources of dust at hi gh red shifts 
(jTodini fc Ferrarall200ll: iNozawa et alll200ll2008[) . 

The origin of the observed dust is not always obvi- 
ous. For SN 2005ip, optical spectra show a progressive 
attenuation of the red wing of both the broad and inter- 
mediate lines, directly confirming the formation of new 
dust in both the ejec ta and post-shocked cool, dense shell 
(jSmith et al.ll2009"bl ). Such direct evidence, however, is 
rare. No more than a handful of supernovae (Lucy et alj 
19911: iMeikle et all 1993: Elm hamd i ct al. 2004, 200l 
Sugerman et al.l 1200 6: Kotafl [200l iPozzo et all [2004 : 



Smith et al.ll2009bn2008b) show direct evidence of dust 



formation i n the ejecta, and a side from SN 2005ip, only 
SNe 1998S (IPozzo et al.l 12004 ') and 2006jc (ISmith et al.l 
[2008bil show direct evidence of dust formation in the cool. 



Table 1 

tzer Observations Summary 



Target Position Flux Int/ Total Ramp X # of 

Field RA Density Bands Pixel Duration Cycles AORS 

DEC (sees) (sees) (sec X #) 



SN 2005ip 09:32:06 0.4 mjy 
+08:26:44 (8/im) 



SN 2005ip 



09:32:06 
+08:26:44 



0.4 mJy 

(8/im) 



IRAC 
all 

IRS 

Spectral Mapping 

low short both 



300 



1062 



9516 



1 



60 X 12 



dense shell (although this scenario is invoked to explain 
the asymm etries observed in the spectr a of SNe 2004 et 
(jKotak et al. 2009) and 2007od (Andre ws et'^ |2010D). 
The observed dust yields, however, all tend to be 2-3 or- 
ders of magnitude lower than required to account for the 
large amounts of dust observed at high redshifts. 

Mid-infrared (3 < A < 15 /xm) observations span the 
peak of the thermal spectral energy distribution from 
dust with temperatures ranging 100 J; Td ^ 1000 K, pro- 
viding strong constraints on the dust mass, temperature, 
and, thereby, the luminosity. These quantities serve as 
useful diagnostics for disentangling the origin and heat- 
ing mechanism of warm dust. Late-time mid-infrared ob- 
servations of supernovae, particularly the Type Iln sub- 
class, are rare. 

In this paper, we present the first late-time (day 936 
post-discovery) Spitzer /1Y{S) spectra of a Type Iln super- 
nova, as well as coincident S pit z er /lYiKG photometry. 
We also present 0.9-2.5 /xm, R = 3000 spectra obtained 
with APO 3.5-m/TripleSpec. Section [2] presents the ob- 
servations and data reduction techniques. The combined 
spectra show evidence for two independent dust com- 
ponents: a hot, near-infrared (HNI) and warm, mid- 
infrared (WMI) component. For each component, we 
derive the dust composition, mass, temperature, and lu- 
minosity. Section|3]explores the origin and heating mech- 
anism of these components to determine the degree to 
which SN 2005ip forms new dust. The HNI dust mass 
originates predominantly from newly formed dust in the 
ejecta, while the WMI component likely originates from 
an 'IR echo' in a pre-existing dust shell. We use these 
results to explore the progenitor system and its evolu- 
tion. Section [5] presents a summary of the findings and 
a discussion of future work. 

2. OBSERVATIONS 

2.1. Spitzer 

As part of PI P 50256, the Spi tzer Infrared Spec- 
trograph (IRS) (jHouck et al.l I2004D obtained one mid- 
infrared spectra on June, 3 2008 (936 days post- 
discovery) with the Short-Low module (SL, R'^60-120, 
5.2-14 /im). The S pitzer Infrared Array Camera (IRAC) 
(jFazio et al.ll2004l ) followed the IRS observations with im- 
ages of the supernova and host galaxy in all four bands 
on June 10, 2008. Table [T] lists the observational details. 

2.1.1. IRAC 

Spitzer collected 5-minute IRAC exposures consist- 
ing of ten 30-second integrations. Pipeline-reduced and 
calibrated (BCD) images were taken from the Spitzer 




Figure 1. Spitzer/TRAC and IRS observations of SN 2005ip in 
NGC 2906 obtained in June, 2008 (~943 and 936 days post- 
discovery, respectively). Showrn here is the Post-BCD 3.5 /im IRAC 
image with the IRS map overlay. SN 2005ip is directly above 
the galaxy, while the galaxy nucleus also falls within the mapping 
scheme. 

archive, and combined into single frames with enhanced 
pixel resolution of 0'.'75 pix"-'^ using the MDPEX software 
package provided by the Spitzer Science Center. Figure 
[T] shows a post-BCD 3.5 ^m IRAC image with the IRS 
map overlay. Since the supernova lies near the center 
of the NGC 2906, the rapidly- varying background of the 
host galaxy complicates aperture photometry. Instead, 
a number of unsaturated, linear, and isolated (e.g. no 
other sources in the wings) stars were used to build an 
empirical PSF for each channel, each of which was used 
to measure the brightness of the supernova. In each chan- 
nel, the residuals from subtracting the best-fit PSF were 
small compared with the predicted uncertainty that the 
photometry task allstar provides, which factors in Pois- 
son noise along with flat-field and profile-fitting errors as 
well as read noise. PSF-fit measurements of field stars in 
the frame were consistent with aperture photometry to 
within 5% in all channels. 

2.1.2. IRS 

The Spitzer /IRS mapped the position of SN 2005ip, 
with 12 cycles of 5 explosures, stepped 2'.'7 perpendicular 
to the 3'.' 7 wide slit. The target spectrum was extracted 
and calibrated from the central pointin g observation 
using the SMART data analysis software (iHigdon et al.l 
120041) ■ with bad pixels identified by IRSCLEAN. Sub- 
tracting off-order observations removed sky and zodia- 
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Figure 2. Spatial profile siiowing SN 2005ip centered in the slit 
with the galaxy nueleus at the left shown by the AdOpt tool within 
the SMART data analysis package. The two point sources are fit 
simultaneously along with the underlying background. 

cal background. All data collection events for a given 
order were then median combined and the spectra ex- 
tracted usin£_y3fi_advfaiiced_optimal extraction routine 
AdOpt (|Lebouteiller et all [20091) . The AdOpt tool pro- 
vides the powerful ability to simultaneously fit multiple 
sources using a super-sampled PSF plus a complex back- 
ground at each row. This feature disentangles the super- 
nova emission from that of the galaxy nucleus as well as 
removes background emission of the galaxy's spiral arms 
(Figure [2]). 

Figure [3] shows the resulting IRS spectrum redshifted 
to account for the radial velocity of the galaxy (2140 
km s~^-|- Vlsr), along with the IRAC photometry. Con- 
tinuum emission peaking at around 3-4 /ini tends to dom- 
inate the spectrum. Few, if any, emission lines are ap- 
parent. 

2.2. TripleSpec 

TripleSpec, an 0.9-2.5 Atm, R = 30 00 spectrograph op- 
erating at APO ([Wilson et al.ll2004t iHerter et al.ll2008[) 
obtained a spectrum on day 862 post-discovery. Forty 
minutes of on-source integration consisted of 8 indepen- 
dent 5 minute exposures nodding between 2 different 
slit positions. We extract the spe ctra with a modifie d 
version of the IDL-based SpexTool (jCushing et al.ll2004 ). 
The underlying galactic arm and sky are approximated 
in SpexTool by a polynomial fit and subtracted from the 
supernova. TripleSpec observations on day 1243 post- 
discovery suggest little evolution occurred between the 
two epochs. We therefore create a single spectrum from 
the near- and mid- infrared spectra (see Figure [3]). 

2.3. Dust Composition, Temperature, and Mass 

Assuming only thermal emission, the combined near- 
and mid-infrared spectra provide a strong constraint on 
the dust mass, temperature, and thereby, the luminos- 
ity. The luminosity of a single spherical dust particle of 
radius, a, and temperature, T^, is given as 



Ld{\) ^ ATia^{nB,{Td)Q,{a)), 



(1) 





>.(^m) 

Figure 3. Combined APO/TripleSpec near-infrared spectra and 
Spitzer/IKAC and IRS mid-infrared data. Although the epochs 
do not match perfectly, TripleSpec observations on day 1243 post- 
discovery (not shown) suggest little evolution occurred since the 
day 862 spectrum plotted here. We therefore treat the combined 
spectra as a single spectrum. The thermal emission that dominates 
the spectrum confirms the presence of warm dust. 



iSmith et al.l (|2009bD provide two pieces of evidence for 
optically thin dust: 1) r « i—^^W ~ 0.5 and 2) 

-^ ■^ ^ ^ ' ^ ^ -'JNIR + -LOPT ^ ' 

the relatively high transmission of X-rays responsible for 
ionizing the unshocked circumstellar medium. For op- 
tically thin dust with mass, Md, at a distance, d, from 
the observer, thermally emitting at a single equilibrium 
temperature, the total fiux can be written as: 



MdB,{Td)K,{a) 



d^ 
where K^{a), the dust mass absorption coefficient, is: 



(2) 



''"^"^"(i^^^)^™'^"^"^^' 



(3) 



where B^{Td) is the Planck blackbody function and 
Qu{a) is the emission efficiency. ,Fox et al., (|2009,) and 



for a dust bulk (volume) density p. 

Given simple dust populations of a single size com- 
posed entirely of either silicate or graphite, FigurelUplots 
the dust absorption coefficient and emission efficiency for 
several grain sizes of each composition, which is derived 
from Mie theory. For the observed spectrum in Figure [3l 
Figure [5] shows the best fit of equation [2] with graphite 
and silicate models using IDE's MPFIT function. The lack 
of an emission feature at ~9 /xm immediately rules out 
any silicate grain contribution. We therefore use only 
graphite models throughout the rest of this paper. 

We fit a multi-component model to the combined spec- 
trum, where each component assumes a single dust mass, 
Md, at a single temperature, Td, composed of graphite 
of a single grain size. The dust mass, temperature, and 
grain size of each component are all free parameters. Fig- 
ure [5] shows the optimal fit, which consists of two domi- 
nant components: a "hot" (^800 K), near-infrared (HNI) 
and "warm" (^^400 K), mid-infrared (WMI) component. 

Adding additional components to the fit tends to yield 
relatively small dust masses that do not improve x^ by 
more than a couple percent. This result has two im- 
plications. First, two unique dust components exist, as 
opposed to a single dust component with a continuous 
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Figure 4. The dust opacity coefBcient, k, and emission efficiency, Q, for both graphite and silicate for several grain sizes as a function of 
wavelength. The emission efficiency is given by the equation [3] The dust opacity coefficient is nearly independent of grain radius assuming 
grains with radii less than the emitting wavelengths (i.e. a < 1 /im at infrared wavelengths). For large grains (a > 1 /im), the dust can be 
approximated by a blackbody (Q = 1). 



temperature distribution. Second, a hotter (> 1000 K) 
third component, if it exists, does not significantly con- 
taminate the near-infrared fits at these late times. (At 
early times, a hot photospheric component is expected 
and is, in fact, ob s erved to dominate the early J-band 
data in iFox et al.l ()2009l ).) We therefore only consider 
two-component fits throughout the rest of this article. 

Table [2] lists the dust masses, temperatures, and lumi- 
nosities associated with various grain size combinations, 
but overall, grain size has little consequence on either 
these parameter values or the resulting goodness-of-fit 
(x^). The only size dependent variable in equation [5] is 
K, but Figure U shows the dust opacity coefficient for 
graphite is independent of grain radius at infrared wave- 
lengths (>1 [iva) for grain sizes <1 /im, which is typical 
for most grains. Other methods (discussed in Sections 
13.21 and l3.3p . however, can constrain the grain size. 

Figure [5] plots the luminosities of both the HNI 
and WMI components on day 936 for 0.1 /im grains, 
along with th e time-series ev olution of SN 2005ip at 
near-infrared (iFox et aLl 120091) and visible luminosities 
(jSmith et all l2009b[) . iSmith et al.l (|2009b[) show that 
as the photospheric component drops off over the first 
~100 days, an optical luminosity plateau, Lpiatoau, ap- 
pears and continues throughout the extent of the obser- 
vations. This plateau arises from radiation generated by 
continuous shock interaction with the dense circumstel- 
lar medium, as opposed to an optical light echo powered 



by the peak supernova luminosity. 

Lacking mid-infrared observations, the data obtained 
prior to day 936 cannot distinguish between the multi- 
ple dust components. The HNI component in Figure [S] 
however, dominates the near- infrared observations (i.e., 
little contribution from the WMI component at near- 
infrared wavelengths) , which plateau throughout the ex- 
tent of the observations (see FigurelS]). In fact, additional 
epochs of TripleSpec spectra obtained on days 862 and 
895 show little evidence for spectral evolution at later 
times. Given these observations, we assume throughout 
this paper that all near-infrared results predominantly 
represent the HNI component evolution. For both the 
HNI and WMI components, the thermal emission arises 
from warm dust. Disentangling the composition, ori- 
gin, and heating mechanism of these dust components, 
however, requires a detailed analysis of potential heating 
mechanisms. 

3. ANALYSIS: DUST ORIGIN AND HEATING MECHANISM 

Section [T] summarized likely origins and heating mech- 
anisms for late-time infrared dust emission, distinguish- 
ing sce narios involv i ng newl y formed versus pre-existing 
dust. ISmith et al.l (|2009bl ) provide spectroscopic evi- 
dence for new dust condensation in the ejecta from days 
60-170 and in the cool, dense shell of post-shocked gas 
at days >413. This dust formation timeline closely cor- 
responds to the evolution of the near-infrared light curve 
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Table 2 

Fitting Parameters 



Fit 



(K) 



That 

(K) 



-ii'-'warm 

(M©) 



Mhot 
(Mq) 



(L©) 



(L©) 



1.0 /im(warni) + 
0.3 /im(warm) + 
0.1 /im(warm) + 
0.01 ^m(warm) - 
0.001 /im(warm) 
1.0 /im(warm) + 
0.5 /im{warm) + 
0.3 /im{warm) + 
0.1 /im(warm) + 
0.01 ^m(warm) - 
0.001 /jm(warm) 
1.0 /im{warm) + 
0.3 /ini(warm) + 
0.1 /im(warni) + 
0.01 ^m(warm) - 
0.001 /ini(warm) 



0.3 /im(hot) 
0.3 /im(hot) 
0.3 /im(hot) 

- 0.3 /im(hot) 
+ 0.3 /.tm(hot) 
0.1 /im(hot) 
0.1 /im(hot) 
0.1 /im(hot) 
0.1 /im(hot) 

-0.1 iim{hot) 
+ 0.1 /.tm(hot) 
0.01 /im{hot) 
0.01 /im(hot) 
0.01 /im(hot) 

- 0.01 /.tm(hot) 
+ 0.01 iim{hot) 



568 
422 
467 
487 
487 
540 
445 
408 
453 
472 
472 
540 
408 
453 
472 
472 



1092 
1061 
1106 
1113 
1113 
847 
836 
838 
859 
862 
862 
897 
887 
910 
913 
913 



5.3 X 
2.6 X 

4.1 X 

4.2 X 
4.2 X 
5.9 X 

1.2 X 

2.8 X 

4.3 X 
4.5 X 
4.5 X 

5.9 X 
2.8 X 
4.3 X 
4.5 X 
4.5 X 



10"'^ 
10-2 
10-2 
10-2 

10-2 

10-^ 

10-2 
10-2 
10-2 
10-2 
10-2 
10-3 
10-2 
10-2 
10-2 
10-2 



4.5 X 

5.7 X 

4.1 X 
3.9 X 
3.9 X 
5.9 X 

6.8 X 
6.7 X 

5.2 X 
5.0 X 
5.0 X 

6.4 X 

7.3 X 
5.7 X 

5.5 X 
5.5 X 



10-5 
10-5 
10-5 
10-5 
10-5 
10-4 
10-"* 
10-4 
10-4 
10-4 
10-4 
10-4 
10-4 
10-4 

10-4 
10-4 



1.2 X 10** 4.2 X lO'^ 
1.1 X 10* 4.6 X lO'^ 

1.3 X 10* 4.1 X lO'^ 
1.3 X 10* 4.0 X lO'^ 
1.3 X 10* 4.0 X lO'^ 
1.0 X 10* 5.8 X 10'^ 
1.0 X 10* 6.2 X lO'^ 

1.0 X 10* 6.2 X 10'' 

1.1 X 10* 5.5 X lO'^ 
1.1 X 10* 5.5 X 10^ 
1.1 X 10* 5.5 X lO'^ 
1.0 X 10* 5.8 X 10^ 

1.0 X 10* 6.2 X 10^ 

1.1 X 10* 5.6 X 10^ 
1.1 X 10* 5.5 X 10^ 
1.1 X 10* 5.5 X 10^ 
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Figure 5. Best fits of tlie grapliite and silicate models given by 
equation [2]to tfie combined near- and mid-infrared data. The lack 
of an emission feature at ~9 fj,m immediately rules out the silicate 
model. The combined spectrum is best fit by a multi-component 
model, consisting of both a "hot" (~800 K), near-infrared (HNI) 
and "warm" (~400 K), mid-infrared (WMI) component. 

in Figure |6l The relative contribution of this newly 
formed dust to the HNI and WMI components, how- 
ever, is not immediately clear, but can be addressed by 
the Spitzer mid-infrared photometry and spectroscopy 
presented here in the context of both pre-existing dust 
scenarios: shock heating and an infrared echo. 

3.1. The Blackbody and Shock Radii 

For both shock heating and an infrared echo, the 
blackbody and shock radii respectively serve as use- 
ful reference points. The blackbody radius, given 

as rbb = ( . "^'14 ) ^ where a is Stefan-Boltzmann's 

constant, defines the minimum shell size of an ob- 
served dust component. In the case of SN 2005ip, 
blackbody fits [Q = 1) of the combined spectrum 
on day 936 in Figure [5] yields blackbody radii of 
r66(WMI) «4.8 X 10^6 cm (0.048 ly) and rb6(HNI) w7.7 x 
10^5 cm (0.0077 ly). 
The shock radius, given as r^ = Vgt for a constant 



JHK-band (1-2.5 urn) < 
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Figure 6. The evolution of SN 2005ip since discovery. Included in 
this plot are near-infrared photometry from Fox et al. (2009), un- 
filtercd (i?-band) photometry from .Smith ct al.. (.2009b.) . and the 
WMI and HNI luminosities on day 936 post-discovery derived from 
the best fits in Section 12.31 and shown in Figure [S] The radioac- 
tive decay for a typical (0.1 M© of 56Ni) supernova is ovcrplotted 
(dashed line). 



velocity, v^, defines the maximum radius that the for- 
ward shock can travel in a time, t. The shock ve- 
locity can be derived fr om the optical emis sion line 
widths. For SN 2005ip, iSmith et all ()2009b[ ) observe 
broad emission lines corresponding to radial velocities 
of ~15,000 km s-i (0.05c) through -900 days post- 
discovery, yielding a maximum shock radius on day 936 
of Tsi ~ 1.25 X 10^'^ cm (0.125 ly). At the same time, 
the intermediate width lines correspond to slower shock 
velocities of -1 000 km s"^ (0.003c) throu gh -900 days 
post-discovery. iChugai fc Danzigeii (ll994j) propose that 
two unique shock velocities can coexist if the progenitor's 
wind is clumpy or asymmetric, as opposed to homoge- 
neous and spherical. A relatively rarefied wind allows 
uninhibited shocks to maintain the observed high veloc- 
ities, while a much slower shock propagates through the 
denser regions. 
Since a shock will typically destroy any pre-existing 
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Figure 7. Post-shock equilibrium dust temperatures (T^j (K)) as a function of post-shock electron density, rio, and temperature, Tc- For 
typical post-shock electron temperatures, Te si lO'^ K, the grains become increasingly transparent to the incident electrons, and the dust 
temperature is only a function of the electron density. 



dust, understanding the shock propagation is essential 
for modehng the evolution of any late-time infrared emis- 
sion. If, for example, the denser regions are concen- 
trated in a homogeneous equatorial disk, only the slower 
shocks need be considered as they propagate contigu- 
ously throughout the disk. Alternatively, if the denser 
regions are concentrated in clumps, the situation is more 
complicated. The fastest shocks will reach the clumps 
first, at which point the shock velocity drops. The aver- 
age shock velocity is a function of the clump filling fac- 
tor and distribution. Since the distribution of the dense 
regions is not well known, we consider a second shock 
radius for which we adopt a more modest shock veloc- 
ity of Vs ~5000 km s~^, which is consistent with the 
inter mediate component of many Type Iln supernovae , 
(e.g. 'Salamanca et al.| 120021: iSmith et al.l l2008al 120101: 
[Steele et al. 2008). This velocity yields an upper limit 
on a second shock radius on day 936 is rs2 ~ 4 x 10^^ cm 
(0.04 ly). 

Both the HNI and WMI shells are optically thin, 
thereby confirming the assumptions presented in Section 
The optical depth of each shell can be written as. 



47^7-2 



(4) 



where Kavg = 435 cm^ g ^ is the absorption coefficient 
for graphite averaged over 1 — 15 /xm. For 0.1 /im grains. 
Table H shows Ma{HNI) « 5.2 x lO"'* Mq and 
Md{WMI) K. 4.3 X 10-2 Mq. Given the minimum 
radii, r > ri,b{HNI) w7.7 x lO^^ cm (0.0077 ly) and 
r > nbiWMI) w4.8 X 10^*^ cm (0.048 ly), equation H 
yields t < 0.6 and 1.3, respectively. 

3.2. Shock Heating 

In the shock heating scenario, hot electrons in the post- 
shock e nvironment collisi onally heat pre-existing dust 
grains. iDwekl (|1987D and iDwek et all ([2008) provide a 
detailed description of this process for silicate dust grains 
and present post-shock equilibrium dust temperatures as 
a function of post-shock electron density, Uc, and tem- 
perature, Tc. Figure [7] presents a similar analysis for 
graphite grains. 

For shock heating to occur, the forward shock must 
have sufficient time to reach the pre-existing dust grains. 



As discussed in further detail in Section [373l the peak su- 
pernova luminosity of SN 2005ip {Lpeak ~ 10^ Lq) vapor- 
izes all dust grains out to a radius, rovap ~ 10^^ cm (0.01 
ly). For the fastest observed shock velocities {vs ~15,000 
km s"^), a shock would require ~70 days to reach the 
evaporation radius, which is consistent with the earliest 
observations of the HNI component (see Figure ^ . The 
earliest observation of the WMI component does not oc- 
cur until the Spitzer observations on day 936, by which 
point even the slower shocks would have crossed the evap- 
oration radius. 

An independent calculation of the shocked dust mass 
provides a useful consistency check. Assuming the hot, 
post-shocked gas heats the dust shell, estimates of the 
total gas mass will determine the dust-to-gas mass ratio. 
The upper limit on the volume of the emitting shell is 
given at 

V;hcn = 47rr2Ar,, (5) 

where Vg is the radius of the shock with velocity, Vs, at 
an age, t, and Ar^ is the shell thickness defined by the 
distance traveled by the shock over the grain sputtering 
hfetime, r^putt, 



An 



1 



r'^s'^sputt: 



(6) 



provided that Tsputt < (^, t-cooi), where TcooI is the radia- 
tive cooling time-scale. The factor of ^ comes from the 
shock jump conditions. For typical pos t -shock gas tem- 
peratures (> 10^ K), IDwek fc Arendtl (|1992[ ) give the 
sputtering lifetime for a grain size, a, and gas density. 



T'sputt(yr) ~ 10'^ 



a(/im) 



(7) 



The total gas mass of the emitting dust shell is therefore 



Afg = ngTOnVshcii- 
Combining equations [5] - [8] yields 



(8) 



Mg(M0) w 8.3 X 10" 



lOOOkms^V \yr) VMm/' 

(9) 
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Table 3 

Shock Heating Mass Predictions at t 
Vs = 15,000 km s-^ 



■ 936 Days for 





WMI 


HNI 


a (/im) 


M, (Mq) 


Md (Mq) 


Md (Mq) 


Md (Mq) 




Predicted 


Observed 


Predicted 


Observed 


0.01 


1.8e-4 


4.5e-2 


1.8e-4 


5.7e-4 


0.1 


1.8e-3 


4.3e-2 


1.8e-3 


5.2e-4 


0.3 


5.5e-3 


2.8e-2 


5.5e-3 


4.1e-5 


0.5 


9.2e-3 


1.2e-2 


- 


- 


1.0 


1.8e-2 


5.9e-3 


- 


- 



which reveals the mass of the shocked gas is independent 
of the grain density. 

Assuming a dust-to-gas mass ratio expected in the 
H-rich envelope of a massive star, Zd = ^ ~ 0.01, 
gives the expected dust mass. For the maximum ob- 
served shock velocity, Vs < 15,000 km s~^, and age, 
t = 936 days. Table [3] compares the predicted dust mass 
to the observed mass listed in Table [2] Assuming shock 
heating is also responsible for the HNI component at 
early times, Table |4] compares the predicted HNI dust 
mass at an age t ~ 70 days to the observed mass listed 
in Table [2j Although no measurement of the dust mass 
exists on day ~70, the relatively const ant HNI lumino s- 
ity (see FigurelH) and temperature (see lFox et al.l (|2OO90 ) 
suggest a relatively constant HNI mass throughout the 
extent of the observations. 

For both the WMI component on day 936 and the HNI 
component on day 70, only large grains (a > 0.3 fj,m) can 
reproduce the observed dust masses. These grain radii 
are large compare d with typical gra in sizes observed in 
supernova shocks (|Dwek et al.l 120081 ). Furthermore, up- 
per limits were assumed for both the sho ck velocity (see 
Secti on [3TT1) and dust-to-gas mass ratio (jWilliams et al.l 
|2006| ) . Lower values would require even larger grain sizes 
to reproduce the observed dust masses. These results 
likely rule out shock heating. 

3.3. Possible Emission From an IR Echo 

For an infrared echo scenario, the supernova luminosity 
heats a shell of dust at a radius, r, to a peak tempera- 
ture, Td- This outer shell may be pre-existing at the 
time of the supernova explosion or it may form when the 
peak supernova luminosity creates a vaporization cavity. 
In either case, light travel time effects cause the ther- 
mal radiation from the dust grains to reach the observer 
over an extended period of time, thereby forming an 'IR 
echo' (Bode & Evans 1980; Dwek 1983). The infrared 
luminosity plateau occurs on year long time scales, cor- 
responding to the light travel time across the inner edge 
of the dust shell. As dust cools from the peak temper- 
ature, it will contribute flux at longer wavelengths. As 
noted in Section 12.31 however, the SN 2005ip spectrum 
is best fit by two components, as opposed to a contin- 
uous temperature distribution. Therefore, this analysis 
assumes a simple light echo model that is dominated by 
flux from only the warmest dust with a single tempera- 
ture, Td. 

The equilibrium dust temperature is set by balancing 
the energy absorbed and emitted by the dust grains. 



Table 4 

Shock Heating Mass Predictions 

at i = 70 Days for Vs = 15,000 

km s~^ 

HNI 
a (Mm) Ma (Mq) M^ (Mq) 
Predicted Observed 



0.01 
0.1 
0.3 



l.Oe-6 
l.Oe-5 
3.1e-5 



5.7e-4 
5.2e-4 
4.1e-5 



where, for a single dust grain. 
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B,{Ts^)Q^Uv)dv (11) 



and 



irad = 47ra^ / 'n:By{Td)Qahs{v)dy 
= —TTpa^ / Bi,{Td)n{v)di' 



(12) 



where rsN is the effective supernova emitting radius, 
Lboi is the UV-optical luminosity, and TgN is the effec- 
tive supernova blackbody temperature, iboi follows from 
Equations [lOl HH and [11 
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par^aT^^ 



j B,{Td)K{v)dv 
J B„{Tsjq)Qe,hs{i^)dv' 



(13) 



^boi depends on the grain radius because although the 
dust opacity coefficient, k, in Figure U is independent of 
grain radius for thermal emission at longer wavelengths, 
it does depend on grain radius for absorption at shorter 
wavelengths (e.g., UV and optical). 

Using equation [131 Figure [8| plots contours of the shell 
size, r, as a function of both luminosity, Lboi, and ob- 
served dust temperature, Td. The luminosity is treated 
as a central point source, assuming the emitting region 
is internal to a spherically symmetric dust shell. Sev- 
eral grain sizes are considered for dust with a graphite 
composition. Although the calculation assumes TgN ~ 
10,000 K, the result is fairly insensitive to this choice. 
The vertical lines show the observed graphite dust tem- 
peratures listed in Table[5]and the approximate vaporiza- 



tion temperature of graphite dust, Te, 



2000 K. The 
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rad7 



(10) 



horizontal lines show both the observed peak, Lpoak, and 
late-time optical/infrared plateau, Lpiatoau, luminosities 
from Figure [6l 

The shaded regions highlight the dust shell radii al- 
lowed by the constraints. The shock radius sets the 
lower limit as any dust within this radius would be in- 
dependently heated or destroyed by the forward shock 
(see Section [321 above). (For the WMI component, the 
blackbody radius, ri,b{WAII), actually sets the minimum 
radius, as described in Section [SHJ) Both the fast, Vgi, 
and slow, rs2, shocks described in Section [3.11 are con- 
sidered, distinguished by the hashed region. The near- 
infrared plateau time-scale sets the upper limit. Figure [6] 
shows the plateau extends for at least ~2.6 years and still 
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Figure 8. Con tour s of the dust shell size, r, plotted as a function of both the supernova luminosity, Lboli and dust temperature, T^, 
given by equation ilS! Several grain sizes are considered for dust writh a graphite composition. The vertical lines shovir the observed graphite 
dust temperatures listed in Table [2] and the approximate vaporization temperature of graphite dust, Tovap ~ 2000 K. The horizontal lines 
show both the observed peak, Lpcak> and late-time optical plateau, ipiatcam luminosities from Figure [B] The shaded regions highlight the 
dust shell radii allowed by the constraints. The shock rad ius sets the lower limit as any dust within this radius would be independently 
heated or destroyed by the forward shock (see Section 13.211 . Both shock radii described in Section |3. H are considered. The horizontal line 
fill distinguishes TsI from rs2- The near-infrared plateau time-scale sets the upper limit. Figure |6] shows the plateau extends for at least 
~2.6 years and still shows little sign of declining. The minimum upper limit is therefore 1.3 light years, although it may certainly be larger. 
Section 13.31 considers three possible echo scenarios. 
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shows Httle sign of dechning. The current upper hmit is 
therefore set by a "plateau radius," Vp, of 1.3 Hght years, 
although the final upper limit will only be determined 
once the light-curve plateau begins to decline. Should 
the plateau extend even longer than 2.6 years the larger 
implied radius only makes a light echo less likely. 

Three possible light echo scenarios exist: 
1) In the first scenario, a sphere of dust completely en- 
compasses the progenitor at the time of the explosion. 
The peak luminosity vaporizes all dust within a radius, 
?'cvap , and warms the inside of the remaining dust shell to 
nearly the v aporization temperature (Tovap ~ 2000 K). 
Figure 8(a) shows that the observed peak optical lumi- 
nosity, -Lpcak, yields a vaporization radius of rovap ~ 
0.013 light years for a = 0.01 ^m graphite grains and 
an e ven s maller radius for larger grains (see Figures |8(b)| 
and 8(c)). These small evaporation radii, however, are 
inconsistent with the observations as a shell of this size 
cannot produce an infrared echo on 3 year time scales. 
Furthermore, these evaporation radii are smaller than 
both shock radii. Even if the actual peak luminosity were 
a factor of 5 larger than observed (Lost '~ 5 x 10 Lq), 
the vaporization radius would be insufficient to account 
for the observations. We therefore rule out this scenario 
for both the HNI and WMI components. 

2) The dust shell inner limit need not lie exactly at 
the vaporization radius. If the progenitor underwent an 
eruption many years before the supernova, the dust shell 
may lie at larger radii. In this second scenario, the peak 
luminosity, Lpeaki heats the dust shell inner radius to 
only the observed temperature. The light echo duration 
therefore defines the minimum cavity radius (i.e., rp—1.3 
light years) . Figure [8] shows that a minimum peak lu- 
minosity of iboi > 5 X 10^" L0 is required to heat a 
dust shell of 1.3 light years in radius to the observed 
WMI temperature, while a minimum peak luminosity of 
Lhoi > 5 X 10^^ L0 is required for the HNI component. A 
larger shell radius requires even larger peak luminosities. 

These required peak luminosities are significantly 
larger than the observed peak luminosity. The true peak 
luminosity of SN 2005ip, however, is not well constrained. 
The earliest R-band photometry was obtained 14-days 
post-discovery, and the discovery occurred a few weeks 
following the actual explosion so that the peak luminosity 
was likely several times larger than the observed early- 
time _R-band luminosity (icst ~ 5 x 10^ L©)- Still, while 
a significant amount of optical absorption might be ex- 
pected by large amounts of p re-existing dust, sig nificant 
reddening was not observed (jSmith et al.l l2009"bl ) and a 
peak luminosity > 10^" Lq is unlikely from this extrap- 
olation. 

The shock breakout in the minutes to hours follow- 
ing the sup ernova explosion may reach peak lumino sities 
> 10" L0 (jSoderberg et al.l[200l Hest et al.l[200l . but 
no such breakout was observed for SN 2005ip. A peak 
luminosity > 10^^ Lq would have made SN 2005ip one 
of the most luminous core-collap se events ever observed 
(iQuimbv e t al. 20071: iRest et al.) '2009'). but figure 11 of 
iRest et all (2009.) suggests SN 2005ip is nearly an order of 
magnitude fainter than SNe 2006gy 2008es, 2005ap, and 
2003ma, especially at early times. Furthermore, opti- 
cal emission from the late-time circumstellar interaction 
successfully accounts for the observed dust temperatures 



(see below). These reasons rule out a light echo driven 
by the peak supernova luminosity for both the HNI and 
WMI components. 

3) A final scenario considers a pre-existing dust shell 
similar to scenario 2 above, but in this case the shell's 
inner radius is located at an intermediate radius between 
the shock and plateau radii (as defined on day 936). The 
late-time optical emission, Lpiatcau, continuously heats 
the dust shell to the observed temperature. This sce- 
nario is not so much a traditional infrared echo as it 
is a reprocessing of the optical emission by the dust. 
Some authors refer t o this as a circumstellar shock echo 
Gerardv et al.l [20021 ).) If the circumstellar interaction 
occurs on a time scale greater than the light travel time 
across the dust shell, the shell radius does not set the 
infrared plateau length. The observed fiux therefore ac- 
counts fo r the enti re she ll. 

Figures 8(a) and 8(b) show that for both a — 0.01 and 
0.1 /xm graphite grams, ipiatoau can heat a dust shell of 
radius r « 0.01 ly to Thni and a shell of radius r w 0. 05 ly 



to TwMi- For a — 0.5 /im graphite grains. Figure^ 

shows that ipiatcau can only heat a dust shell of radius 
r < 0.005 ly to Trni and a shell of radius r < 0.04 ly to 
TwMi- In this scenario, both the fast {rsi = 0.125 ly) 
and slow (rs2 = 0.042 ly) shock radii are larger than the 
HNI shell radii by day 936, ruling out a light echo of this 
sort for the HNI component. The same is true for the 
WMI shell composed of larger grains (a > 0.5 /zm). In 
the case of the WMI shell composed of smaller grains 
(a < 0.1 /zm), however, the slower shock has not yet 
reached the shell radius (r fv 0.05 ly). Not only is this 
scenario possible for the WMI component, but the radius 
is consistent with the WMI blackbody radius (rbb (WMI) 
= 0.048 ly). Furthermore, the grain sizes are typical of 
those o bserved in other s upernova circumstellar environ- 
ments (jDweket al.ll2008| ). 

Predictions for the light curve evolution can be made 
from this model. Dust that remains at radii beyond the 
slower shock radius, rs2, will continue to radiate and 
contribute to the light echo plateau. Although the ex- 
act distribution of these dense, dusty regions is not well- 
known, this scenario suggests they must be distributed in 
such a way that the fastest shocks do not interact with 
the dust. As described in Section 13.11 these dense re- 
gions may be concentrated in either an equatorial disk 
or clumps. The clumps must have a large filling factor 
if the fastest shocks are not to interact with a significant 
portion of the dust. As the slower shocks continues to 
expand, however, they will ultimately destroy the dust. 
Assuming the dust lies at a radius consistent with 0.1 
^m grains (see parameters in Table [2]) , the WMI flux 

will begin to decrease at i « — ^ 2775 days post- 
discovery and continue to decrease as a function of the 
emitting shock radius and the dust distribution. While 
this scenario is entirely consistent with the WMI obser- 
vations at the current time, mid-infrared observations at 
later epochs can reveal the accuracy of this model's pre- 
dictions. 

4. DISCUSSION 

4.1. The HNI Component 

The above HNI component analysis rules out the pos- 
sibility of both pre-existing dust scenarios (i.e., shock 
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heating or infrared echo). Only condensation models 
(i.e., in the ejecta or cool, dense shell) remain as viable 
scenar ios for expla i ning th e origin of the HNI compo- 
nent. iSmith et al.l (|2009bl ) confirm new dust formation 
via extinction in the broad (^-^15,000 km s~^) Ha wings 
between days ^60-170 and in the intermediate (~2000 
km s^^) He I wings at days >413, but the location of 
this new dust remains ambiguous. The extinction in the 
broad and intermediate components suggest new dust in 
the fast ejecta at early times and post-shock cool, dense 
shell at later times, respectively. 

To form dust at the high ejecta velocities, however, 
is difficult due to the lack of heavy metals traveling at 
these speeds. If the dust formed in the cool, dense shell 
at early ti mes, as suggested in the case of SN 2006jc 
(iMattila et al. 2008), this would ex plain the att e nuation 
of the broad Ha line observed bv ISmith et al.l (|2009bf ) 
without having to invoke dust formation at high veloci- 
ties, but would fail to explain the lack of observed inter- 
mediate wi dth He I line s produced by the circumstellar 
interaction. ISmith et al . ( 2009b) propose one alternative 
scenario in which the dust may form in the post-shock gas 
of individual clumps, which then become incorporated 
into the expanding fast ejecta when a clump is eventu- 
ally destroyed. This scenario, however, requires several 
assumptions, most of which rely on unknown clumping 
properties. While the available observations limit the 
ability to isolate the precise region of dust formation, 
the relatively fiat near-infrared flux in Figure [6] suggests 
a majority of the dust contributing to the HNI flux must 
have formed at early epochs. 

Newly formed dust cannot reproduce the observed HNI 
luminosity plateau solely by cooling from condensation 
to the observed temperature Thni ~ 800 K. Given an 
average energy per particle e = CnAr^ for the specific 



heat for graphite Cg (given bv lDraine fc Lil (|2001f )) and 
AT K, 1200 K, an unsustainable mass condensation rate 
of M = Thni/c ~ 100 M0 day^^ would be necessary 
to reproduce the observed fiux (see a more detailed ex- 
planation of this argument for the case of SN 2006jc in 
iFoxet al.l(l2009l) l. 

Instead, an alternative heating mechanism must power 
the thermal emission from this newly formed dust. Fig- 
ure [6] shows that radioactive heating is insufficient to 
power the late-time near-infrared emission. While heat- 
ing by the reverse shock is possible, a similar analysis as 
performed in Section r3.2l suggests this scenario is unlikely 
because only large grains (a > 0.5 /zm) can reproduce the 
observed flux. More likely, the optical luminosity gen- 
erated from the forward shock interaction continuously 
heats the newly condensed dust in the same way that 
it heats the WMI component discussed in Section 13.31 
Unlike the infrared echo scenario, however, the newly 
formed dust exists interior to the forward shock radius 
and the shock emission cannot be treated as a central 
point source. 

4.2. A Multi- Component Model 

A multi-component dust model for SN 2005ip now be- 
gins to emerge, composed of both an inner, "hot" (^-^800 
K), near-infrared (HNI) and an outer, "warm" ('^400 K), 
mid- infrared (WMI) component. Newly formed dust in 
the either the ejecta or cool, dense shell likely dominates 
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Figure 9. Illustration of the proposed multi-component model for 
SN 2005ip consisting of both a hot, inner (HNI) component and 
warm, outer (WMI) component. The HNI dust mass originates pri- 
marily from newly formed dust in the ejecta, while the WMI com- 
ponent likely originates from an circumstellar shock echo that forms 
from the heating of a large, pre-existing dust shell. CoUisional radi- 
ation from continuous shock interaction with the dense circumstel- 
lar me dium generates the optical luminosity plateau (Smith et al] 
I2009bl ). which is the likely heating mechanism for both the HNI 
and WMI components. 



the HNI component. The WMI temperature and black- 
body radius, r66(WMI) «4.8 x lO^^ cm (0.048 ly), are 
consistent with a pre-existing dust shell heated by the 
observed late-time optical luminosity generated by the 
forward shock interaction (Section I3.3p . Figure JH] illus- 
trates the locations of each component, as well as the 
likely origins and heating mechanisms. 

The large pre-existing dust mass that contributes to 
the WMI flux component suggests significant mass loss 
from the progenitor. Assuming a dust-to-gas ratio 
Zd ~ 0.01, the observed WMI dust mass listed in Ta- 
ble H (Md(WMI) ~ 0.01 - 0.05 Mq) yields a total gas 
mass of Mg(WMI) ^1-5 M©. This mass accounts 
for the entire WMI emission, as the entire shell con- 
tributes to the observed flux given the size of the shell 
(r6h(WMI) «4.8 X 10^*^ cm (0.048 ly)) is significantly less 
than the observational time-scale (936 days). The asso- 
ciated mass loss rate is 



^^M,(WMI) 



Ar 
:7.5x 10" 



.3/ Mg(WMI) x 
V Mq / 

t)(^ 



(14) 
(15) 



{i-> 



yr 



"). 



V 120 km s- 

On day 413 post-discovery, ISmith et al.l ()2009bl ) mea- 
sured a progenitor wind velocity Vyj = 120 km s~^. As- 
suming a constant wind velocity and thin shell (Ar = 
^), the mass loss rate is M « 7.5 x 10^^ - 3.8 x 10"^ 
Mq yr~^, which is two-three orders of m agnitude larger 
than calculated by ISmith et al.l ()2009bf) at inner radii. 
Furthermore, this mass loss rate is likely a lower limit. 
He I P-Cygni profiles at 1.083 jjluy from the TripleSpec 
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Apparent Doppler Velocity (10* km s'') 

Figure 10. The 1.083 /^m He I line. The shocked gas creates 
the broad (~ 1000 km s~^) emission feature. The unshocked, 
slow moving, circuinstellar medium formed by the progenitor wind 
produces the superimposed P Cygni profile. The blue edge of the P 
Cygni absorption feature, ~200 km s~^ , measures the wind velocity 
on day 862 post-discovery. 



spectra show the progenitor wind velocity may have been 
closer to 200 km s~^ at the WMI radius (Figure ITU)) . as- 
suming the P-Cygni feature is generated by the cool, low 
velocity circumstellar environment that coincides with 
the WMI shell. 

The larger mass loss rate (and possibly faster pro- 
genitor wind velocity) of the WMI shell suggests a 
denser progenitor wind (or "eruption") occurred at i < 
!ii(}^ « Y2^l°^' j^f ^ 125 years pr ior to the core 
collapse. Although iSmith et al.l (|2009bl ) conclude that 
the SN 2005ip progenitor was likely a red supergiant, 
these stars typically only have wind speeds •(;„; ~ 20 — 
40 km s~^ and mass loss rates up to M = 10~* — 10"'^ 
Mq yr-i (jSmith et al.ll2009aD . The observed charac- 
teristics associated with the WMI shell are more con- 
sistent with Luminous Blue Variable (LBV) sta rs (e.g., 
IDavidsonI 119891: IHumphrevs fc DavidsoiJ I1994D . which 
can have wind speeds on order of hundreds km s~^ (e.g., 
ILeithereil[T997; Kotak fc Vink 2006) and can have mass 
loss rates up to M = 10"^ Mr:, vr"^ (ISmith fc Owockil 
[200l ISmith fc Hartiganl[200l ISmith et al.ll2007n . 

5. CONCLUSION 

The Spitzer spectrum presented here confirms the 
presence of warm dust in SN 2005ip. Combined with 
near-infrared observations, the results show evidence for 
two independent dust masses: a hot, near- infrared (HNI) 
and warm, mid- infrared (WMI) component. Infrared ob- 
servations span the peak of the thermal emission, thereby 
providing strong constraints on the dust mass, tempera- 
ture, and luminosity, which serve as critical diagnostics 
for disentangling the origin and heating mechanism of 
each component. The HNI dust mass originates primar- 
ily from newly formed dust in the ejecta, or possibly the 
cool, dense shell, while the WMI component likely orig- 
inates from an circumstellar shock echo that forms from 
the heating of a large, pre-existing dust shell. For both 
components, the heating mechanism is likely the opti- 
cal luminosity generated from the forward shock interac- 



tion with the circumstellar medium. For wind speeds of 
- 100 - 200 km s-i, the WMI dust shell likely formed 
via an eruption ~ 100 years before the supernova. These 
characteristics are consistent with a LBV progenitor, 
which has been linked to some core-collapse supernovae 
(Kotak fc Vink 2006 ). but is an emerging trend p articu- 
larly within the Type Iln subclass (e.g .,TGal-Yam et al. 



20071: 'Smith et al."2007l: lTrundle et al]|200& .Smith et al. 
2008a, 2010; Gal- Yam fc Leonardll2009ir 

The newly formed HNI dust mass, M<i(HNI) ~ 5 x 
10""* Mq, is at least two orders of magnitude lower than 
predicted by ejecta condensation models to reproduce 
the large amounts of dust observed at high reds hifts 
(jTodini fc Ferraral [200ll: INozawa et all [20031 l200l . A 
dust mass of ~ 10~^ M(n is comparable to other Type 
Iln events (eg.. IGerardv et al. 2002; Po zzo et al.ll2004l: 
IMeikle et al.l 120071 ). For clumpy eject a, it should be 
noted that both semi-analytical models (jVarosi fc Dwekl 
1 19991) and Monte C arlo radiative-transfer simulations 
( Ercolano et al.|[2007D have shown that dust masses can 
be at least an order of magnitude larger than that pre- 
dicted by the analysis presented in Section 12.31 The 
clump structure for any supernova, however, remains un- 
constrained at present. The pre-existing WMI dust mass, 
Md(WMI) ~ 0.05 Mq, is quite significant. Although, 
in the case of SN 2005ip, the forward shock may ulti- 
mately destroy this dust (mid-infrared observations at 
later epochs will reveal the accuracy of this model), the 
forward shocks may be sufficiently decelerated in other 
systems to allow for dust survival. 

Late-time mid-infrared observations of dust in the su- 
pernova environment provide for a unique interpretation 
of the circumstellar environment and progenitor system. 
Nonetheless, mid-infrared observations are quite rare. 
Presented in this paper is the first mid-infrared spectrum 
of any Type Iln supernova. In the future, we hope to 
grow the database of mid-infrared observations of Type 
Iln supernovae. Doing so, however, can be slow as Type 
Iln events are particularly rare, consisting of only ^ 2-3% 
of all core-collapse supernovae (jGal-Yam et al.ll2007l ) and 
occurring at a rate of no m ore than 10/yr out to 150 Mpc 
(jDahlen fc Franssonlll999l ). We therefore plan to revisit 
all Type Iln supernovae from the past ten years to deter- 
mine the degree to which this subclass exhibits late-time 
dust emission and identify the emission mechanism. 
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